Transition from sub-Poissonian to super-Poissonian shot noise in resonant quantum 

wells 
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We investigate the transition from sub-Poissonian to super-Poissonian values of the zero- 
temperature shot noise power of a resonant double barrier of macroscopic cross-section. This transi- 
tion occurs for driving voltages which are sufficiently strong to bring the system near an instability 
threshold. It is shown that interactions in combination with the energy dependence of the tunneling 
rates dramatically affect the noise level in such a system. Interaction-induced fluctuations of the 
band bottom of the well contribute to the noise and lead to a new energy in the Fano factor. They 
can enhance the noise to super-Poissonian values in a voltage range preceding the instability thresh- 
old of the system. At low voltages, interactions may either enhance or suppress the noise compared 
to the the non-interacting case. 

PACS numbers: 73.50.Td, 73.23.-b, 73.20.Dx, 73.61.-r 



I. INTRODUCTION 



The investigation of the current noise properties pro- 
vides an important tool to gain additional information on 
electric systems. A standard measure of the noise is the 
Fano factor F which is denned as the ratio of the current 
noise spectral density Si (to) and the full Poisson noise 
2e(I) where (I) is the average current driven through the 
well. Thus the Fano factor F — Si(ui)/2e(I) is smaller 
than one for sub-Poissonian noise and exceeds one for 
super-Poissonian noise. In this work we present theoreti- 
cal results for the Fano factor of a resonant double barrier 
with the macroscopic cross-section and investigate the 
transition from su6-Poissonian noise to swper-Poissonian 
noise. 

An early experiment of Li et afJ focused on the sup- 
pression of the noise power below the Poisson value. Al- 
ready in this work it was noticed that the Fano factor 
increases from the suppressed values in the voltage range 
of positive differential conductance (PDC) as the range of 
negative differential conductance (NDC) is approached. 
Experiments by BrownEI and Iannaccone et aa showed 
that in the NDC range, the current noise increases even 
above the Poisson value, and is thus described by * Fano 
factor exceeding one. Recently, Kuznetsov et aa have 
presented a detailed investigation of the noise oscillations 
from sub-Poissonian to super-Poissonian values of a reso- 
nant quantum well in the presence of a parallel magnetic 
field. The magnetic field leads to multiple voltage ranges 
of NDC and permits a clear demonstration of the effect. 
Below we present a theory which describes the transition 
from sub-Poissonian (F < 1) to super-Poissonian noise 

(F>1) ' 

Our approach treats the fluctuations of the potential 
and the charge in the well. The theory describes the 
transition from sub-Poissonian noise to super-Poissonian 
noise at large voltages. At smaller voltages it gives ei- 
ther an additional suppression of the noise power or an 



enhancement of the noise power compared to theories 
which neglect the potential fluctuations. We find that a 
new quantity, an energy A, enters the Fano factor. The 
energy A is a measure of the sensitivity of the noise to a 
collective charge and current response generated by the 
Coulomb interaction. The interaction energy A is deter- 
mined by the response of the average current (/) to a 
variation of the electrostatic potential U in the well 



J = 



du 



and by the response of the charge 



Co 



d(Q) 
dU 



(1) 



(2) 



to the potential in the well. Co has the dimension of ca- 
pacitance. It is positive in the PDC range, giving rise 
to screening of excess charges, and becomes negative as 
the NDC range is approached, corresponding to over- 
screening of excess charges. For a double barrier with a 
geometrical capacitance Cl across the left barrier and a 
geometrical capacitance Cr across the right barrier we 
find an interaction energy 



A = 



hJ 



Cl + Cr + Co 



(3) 



Note that J has the dimension of a conductance. The 
interaction energy thus has the form of a dimensionless 
conductance hJ/e 2 multiplied by an effective charging 
energy of the well e 2 /{C L + Cr + C ). In the NDC range 
the charging energy grows and leads to a very large in- 
teraction energy A. For a resonant level characterized 
by decay widths Tl,r through the left and right barrier, 
an asymmetry Ar = (Tl — Tr)/2, and a total width 
r = + Tr , we find a Fano factor given by 



2 + T 2 ■ 



(4) 
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Note that the escape rates T^r are positive quantities 
but A and Ar can have either sign. For a symmetric 
double barrier the sign of Ar is independent of the ap- 
plied voltage. However, as will be shown in this work, 
for any barrier A changes sign as a function of voltage: 
It is negative in the stable part of the I-V-characteristic 
and becomes positive and large in the NDC part of the 
I-V characteristic. 

If A is set to zero in Eq. (Q), the Fano factor for 
non-interacting electrons of Chen and Tingd is recovered. 
This result is a consequence of partition noise: an elec- 
tron scattered at an obstacle with transmission probabil- 
ity T can only|-.b€ either transmitted or reflected. The 
partition noiseETEl gives rise to shot noise proportional 
to T(l — T). For a resonant tunneling barriena, using 
the Breit-Wigner expression for the transmission prob- 
ability, this leads immediately to the Fano factor, Eq. 
(||) with A = 0. For A = the Fano factor is mini- 
mal (equal to F = 1/2) if the two barriers lead to equal 
escape rates Ar = and the Fano factor is maximal 
(with a value of F = 1), if the double barrier structure is 
very asymmetric, |Ar| ~ T/2. The lower limit F = 1/2 
is obtained, independent of whether transmission is a 
quantum mechanical, fully coherent process cp-flj-sequen- 
tial two step process without phase coherencellHilJ. This 
maximal suppression is also obtained in the presence of 
a magnetic field (see Ref. [t5|). On the other hand, mod- 
els with internal degrees of freedomtJ can give rise to a 
noise suppression slightly below F = 1/2, if the internal 
degrees of freedom couple also to the reservoirs. 

Our result, Eq. (Jl|), also has a Fano factor of 1/2 as 
a lower limit. However, depending on the asymmetry, at 
low voltages our theory predicts either a reduction or en- 
hancement of the Fano factor compared to a free-electron 
theory. A negative A arises if an increase of the internal 
voltage U decreases the current through the well. This 
is the case in the stable range of the I-V characteristic: a 
charging of the well is counteracted by the Coulomb inter- 
action which creates a positive potential response which 
in turn drives the resonant level further up in energy. 
This leads to a Coulomb induced correction of the Fano 
factor for a resonant double well even at low voltages. 

As the range of the NDC is approached, the interac- 
tion energy A changes sign and becomes positive to such 
a degree that it completely dominates the Fano factor. 
This large increase is again due to the Coulomb interac- 
tion. However, now we are in a range where the average 
charge increases in response to a positive variation of 
the internal voltage. A charge injected into the well cre- 
ates a fluctuation of the internal potential which causes 
even more charge to flow into the well. Thus the current 
fluctuations are magnified and give rise to a Fano factor 
which exceeds the Poisson value. 

Large fluctuations are commonly expected near equi- 
librium or non-equilibrium phase transitions. The anal- 
ogy of the far from equilibrium current instabilities with 
equilibrium phase transitions was the subject of work 
by Pytte and ThomasE-J These authors investigated the 



critical slowing down of dielectric relaxation modes at a 
bulk current instability (CLunn effect). Subsequent work 
by Biittiker and Thomasta analyzed a one band-current 
instability (superlattice) and discussed in addition to the 
dielectric relaxation modes also the modes which arise 
due to the coupling of the current fluctuations to the 
magnetic field. In this work we investigate the resonant 
well under the condition of complete external voltage con- 
trol and take into account the electrostatic but not the 
magnetic part of the Coulomb interaction. 

There are two theoretical worksiEl which also indicate 
Fano factors in excess of the Poisson value. However, 
these works do not include the partition noise. Instead, 
it is assumed that the incident current on the cathode 
side is fully Poissonian but that the transmitted and re- 
flected carrier streams are modified through the current 
dependence of the transmission probabilities. Our the- 
ory takes as the basic noise source the partition noise and 
treats the associated charge fluctuations self-consistcntly. 
To treat the potential fluctuations, we evaluate the off- 
diagonal elements of a potential operator of the well. In 
contrast, Ref. ||does not treat potential fluctuations and 
Ref. [l9] includes self-consistant-.effects only on the level 
of probabilities. These worksBtJ also predict in the PDC 
range Fano factors lower than F = 1/2. Such low Fano 
factors have been observed in the recent experiment by 
Kuznetsov et aa. However, since in Refs. pip)] an es- 
sential source of noise, the partition noise, which gives 
raise to the minimal Fano factor F = 1/2 in the non- 
interacting case, is absent, it is perhaps not surprising 
that lower noise limits are predicted than by our theory. 

Below we present a self-consistent, analytical calcula- 
tion of transport properties of resonant tunneling quan- 
tum wells in the limit when the total decay width of 
the resonant state is lower than all other energy scales 
of the problem. We consider the zero-temperature case 
and use the system of units with h = 1. We continue to 
use the expression " shot noise" even in the case when the 
noise is dominated by interaction. We emphasize, that in 
this case, we deal not with single independent electrons 
traversing the conductor, but with a collective response 
of many electrons. The term "shot noise" underlines here 
that the basic source of noise is the partition noise, inde- 
pendent of the degree of interaction. 

The paper is organized as follows: In Section II we in- 
troduce the self-consistency equation and solve it for the 
average quantities. Sec. Ill is devoted to the average cur- 
rent. The results of these two sections by themselves are 
not new, and similar discussion can be found in the liter- 
ature (see e.g. Ref. |l|). However, they are necessary for 
the subsequent calculation of the current-current fluctu- 
ations. This calculation is performed in Sec. IV. In Sec. 
V we present the conclusions and discuss in more detail 
the conditions of validity of our theory. 
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II. AVERAGE CHARGE AND POTENTIAL 



A. The charge of the well and the self-consistency 
equation 

We consider a quantum welfi extended in the x and 
y-directions (to be denoted as _!_) with the area A. The 
potential profile in the ^-direction is shown in Fig. 1. 
We assume that the longitudinal and the transverse mo- 
tion of the electrons are separable. The motion in the z- 
direction is quantized, and there exists a set of resonant 
levels with energies E n , n = 0, 1, . . .. The transverse mo- 
tion is described by a continuum of momenta p± . Thus 
the expression for the total energy (measured from the 
band bottom of the well) is 



R — R 



■p\/2m. 



The density of states in the well can be factorized into 
longitudinal and transverse parts, 

v(E) = v 2 v z {E z ), E = E z +pl/2m. 

Here v 2 = mjl-K is the the two-dimensional density of 
states (per spin). 



R 



u 



Ep 



FIG. 1. The potential profile of the quantum well. 

For interacting electrons the band bottom of the well 
u = eU (see Fig. ([!])) is not an independent quantity but 
must be found from a self-consistency equation. In the 
present work, we account for the interaction by including 
the energy required to charge the well. We introduce the 
capacitances Cl and C R of the left and the right bar- 
rier, respectively. The charge allowed in the well by the 
Coulomb interaction is then determined by these capac- 
itances and the voltage drop across the left barrier u — v 
and across the right barrier u. Here v > is the voltage 
applied across the double barrier structure. (We assume 
that the Fermi energies of the left and right contact are 
equal). Thus the total charge Q[u] in the well is 



Note that the charge Q is always non-negative, and 
thus the band bottom either assumes the value u e = 
(Cl + Cr)~ 1 Clv of the empty well, or lies above it. 

Throughout this work, we describe scattering at the 
double barrier in the Breit-Wigner limit. In particular, 
the longitudinal density of states is 



r(E z 



{E z -E n y+T*{E Z )IA- 



(6) 



Here we have introduced the total width of a single res- 
onance T(E Z ) = T L (E Z ) + T R (E Z ) with T L (E Z ) and 
Fr(E z ) being the partial widths for decay through the 
left and right barrier, respectively. The Breit-Wigner 
limit is appropriate in the limit of small tunneling prob- 
abilities through the two barriers and provides an excel- 
lent description as long as the resonant level is not in the 
immediate vicinity of the band bottom of any reservoir, 
u + E n = v. 

In the standard application of the Breit-Wigner ex- 
pression the partial width is independent of energy; for 
energies E z close to the resonant state with energy E n 
the width is evaluated at the energy of the resonant state, 
T(E n ). Here, in order to describe the density of states 
in the limit where the resonant level approaches the con- 
duction band bottom, we consider the decay widths to 
be energy dependent. We determine the energy depen- 
dence by considering the quantum-mechanical problem 
of transmission through a single rectangular barrier. We 
write 

Tl{E z ) = a L El' 2 (E z +u- vf/ 2 e{E z )e{E z +u-v); 
T R (E Z ) = a B E x J 2 (E z + u) l l 2 d(E z )B(E z + u), (7) 

where and an are dimensionless parameters. 

With the density of states as specified above, we can 
now find the charge Q[u] inside the well. For this pur- 
pose, we divide the total longitudinal density of states 
v z (E) given by Eq. (||) into the partial densities of states 
of electrons which are incident from the left or the right 
reservoirs, 



We fine 



v z (E z )=u L {E z ) + u R (E z ). 



vl,r\E z ) = „,„ » v z (E z ). 



(8) 



With the distribution functions Jl,r of the left and right 
reservoir and the short-hand notation E± = p 2 L /2m J we 
find for the charge of the well 



= e E [ dE * E "i(E z )fi(E z + E 



(9) 



i=L,R 



Cl{u — v) + C R u — eQ[u] 



(5) 



Eq. (^) gives the total charge in the well in terms of the 
carriers injected from each reservoir. 

To proceed we now specify two additional conditions: 
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(i) Of all longitudinal energy levels only Eg is relevant, 

all others lie too high. In addition, we assume for 
definiteness E > Ep. 

(ii) The decay widths T l and T R are much smaller then 

all other energy scales. In particular, we will re- 



place (£ 2 
ate. 



T 2 /4) 1 by 27iT~ 1 <!>(£), where appropri- 



Making use of these conditions and replacing the sum 
over p± by an integral, we obtain 



Q[ u ] = { Ce T T(Eo) ^ + Ef ~ u - E ^> v-E a <u<v-E + E F 



(10) 



0. 



otherwise 



Here we have introduced the "quantum capacitance" 
C e = ev-iA of the two-dimensional electron gas. The 
dependence of Q [u] on u is shown schematically in Fig. 2 
(solid line). 



B. Self-consistent calculation of the band bottom u 

Equation (^) expresses the band bottom in the well u 
as a function of the applied voltage v. The equation is 
cubic and can be solved analytically. However, the re- 
sulting expressions are cumbersome and not very trans- 
parent. Therefore, we discuss now the behavior of the 
function u(v) qualitatively, and support the discussion 
by numerical results (Fig. 2). 

The left-hand side of Eq. (||) turns to zero (empty 
well) for u = (Cl + C r )~ 1 Clv. The charge is non-zero 
only between u = v — Eo and u = v — Eq + Ep. With 
the help of the voltages 

_ c L + c R _ c L + c R 

v a = — [E - E F ) ; v b = E , 

<^R <^R 

we can distinguish the following four regimes: 

I. < v < v a (Fig. 2a). The well is uncharged; the only 
solution of Eq. (g) is, 

C L 



u(v) 



Cl + Cr 



v. 



(11) 



II. v a < v < Vb (Fig. 2b). The well is now charged, but 

Eq. (||) has still only one solution. 

III. Vb < v < v* (Fig. 2c). For voltages above Vb three 

solutions emerge; now ui(v) is given by Eq. (|TT| ) 
and corresponds to the uncharged well, while two 
others U2(v) and us(v) describe a charged well. The 
lesser of them, ii2(v), is unstable, and the greater 
one, us(v), is stable. At a certain voltage v* , which 
is a complicated function of the parameters of the 
system, the straight line e~ 1 {CL{u — v) +Cru] is a 
tangent to Q[u\. This voltage determines the upper 
margin for this regime. 

IV. v* < v (Fig. 2d). The well is again uncharged, and 

the only solution is given by Eq. (p]). 

The dependence u[v] is illustrated in Fig. 3. 



III. AVERAGE CURRENT 



To find the current across the quantum well, we start 
from the general expression for the current density, 



2ttA 



[ dE z T(E z ) 
p± J 



x [f L (E z + E ± ) - f R (E z + E±)] , (12) 
where T(E Z ) is the transmission coefficient, equal toil 
T L (E Z )T R (E Z ) 



T{E Z 



(E z ~Eo) 2 +T2(E z )/4- 



(13) 



Treating Eq. ( |l2"| ) in the same way as previously Eq. (|^), 
we obtain 

C e T L {E )T R (E ) 
3= A—^) iv + E F -u-Eo) (14) 

for v < u + Eq < v + Ep and zero otherwise. 

This can be easily translated into the current-voltage 
characteristics, given the known dependence u(v). The 
current is zero for v < v a and v > v*. For v a < v < Vb 
it is given by Eq. (p~4|) , where u is replaced by the only 
solution u(v). Finally, for Vb < v < v* the three possible 
solutions u(v) yield three values for the current for any 
voltage. One of them is j = and corresponds to the 
uncharged well (the solution ui(v)). Two other branches 
correspond to the charged well; the stable branch (3 in 
Fig. 3, due to 113(1))) lies above the unstable one (2). The 
current-voltage characteristics is shown in Fig. 4 (solid 
line). The differential resistance jumps from zero to a 
finite positive value at v = v a and then decreases with 
voltage, passing through zero at a certain point. It is 
negative for voltages close to v* and turns to —00 at this 
point. 

In the experiment the system is expected to exhibit 
a hysteretic behavior (cf. Ref. [l2]). Upon increasing 
the voltage, the current follows the upper branch (3 in 
Fig. 4) until v = v* , and then jumps onto the zero-current 
branch. Upon decreasing the voltage, the current is zero 
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until v = Vb, and then jumps to a finite value on the 
branch 3. The unstable branch 2 is not observed. Note 
that the range of hysteresis v* — Vb is broader than the 
range of negative differential resistance. 
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(a) 


Q 














(b) 



u 





u 



(d) 



FIG. 2. Solutions of the self-consistency equation (|^). The 
solid line represents the charge Q[u]; the dashed line is 
e~ [Cl{u — v) + Cru]. The position of the vertical axis is ar- 
bitrary. The following set of parameters is chosen: cll = clr; 
Cl = Cr = C e /10; E Q = 3E F /2. Voltages are (in units of 
E ) 1/3 (a), 1 (b), 7/3 (c), 3 (d). The voltage v* is equal 
to 2.70. The solutions are shown as Ui\ in the case (c) there 
are three solutions, of which ui and uz are stable, and U2 is 
unstable. 

The above analysis is valid in the case of low trans- 
parency, T(Eq) -C v* — Vb- In the opposite case the hys- 
teretic range is smeared, and a broad range of negative 
differential conductance appears instead. For weak inter- 
action C e <^Cl, Cr the width of the hysteretic range is 
rather small, 



Vb 



( a L 



XIclrCl+Cr) Eo 



but it increases with the capacitance C e and eventually 
must exceed the smearing T(Eq) of the resonant level. 



u 












3/7 
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IV 



FIG. 3. Dependence u(v) (solid line). Curves 1, 2, and 3 
in the regime III correspond to the solutions Ui, 112, and 1/3, 
respectively. The solution U2 is unstable. The parameter set 
is the same as for Fig. 2. 



IV. FLUCTUATIONS 

A. Charge operator 

The fluctuations of the charge, the potential and the 
current are determined by the off-diagonal elements of 
these quantities. We must use operator expressions in- 
stead of average quantities. Following Ref. ^3|, we intro- 
duce creation and annihilation operators aJ^iS^pjJ and 
<ia(E z ,p±) of carriers in the reservoir a = L,R. The 
quantum statistical average of bilinear products of these 
operators is 



(al(E z ,p ± )a (E' z ,p' ± )) = 5 a p8 p±pl± 
xS(E z -E' z )f a (E z +pl/2m). 



(15) 
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Here we are interested in the zero-frequency spectral 
correlations. In this limit the off-diagonal elements of 
the charge operator can be obtained with the help of off- 
diagonal density of states elements which in turn can be 
related to energy derivatives of the scattering matrix. In 
the presence of a current amplitude incident at contact 
a and a current amplitude at contact (3 the off-diagonal 
density of states .element is a p , in terms of scattering ma- 
trices s a p, isE-JEj 



7 L 



ds 



7/3 



dsl 



7 a 



7Q dE r 



dE, 



(16) 



Using the explicit expressions for the scattering 
matricesEJ, 



s Q/3 = exp [i (4> a + 4>jj)} \d a - i 



E z -E + »r/2 



(17) 



(with cj) a being arbitrary phases), and having in mind 
that the main singularity is generated by the derivative 
of the denominator in Eq. (|17|), we easily obtain 



v a p{E z ) 



1 y/T a (E z )T p (E z ) 
2n {E z -E Q f+Y\E z )/A' 



(18) 



With the help of these off-diagonal density of states ele- 
ments the operator of the charge of the well iscfl 



= e EE / duje ~ Wt [ dE z u a0 (E z ) 

a/3 p± 

x d^ a {E Zl p 1 _)dp{E z + u>,px)- 



(19) 



It is straightforward to check that Eq. (|l^) reproduces 
the average charge (j|). Below we will use the charge 
operator to find the fluctuations. 




FIG. 4. The current-voltage characteristics j(v) (solid 
line). The hysteretic behavior is shown by dashed lines. The 
parameter set is the same as for Fig. 2. 



B. Band bottom operator 

With Q being an operator, the self-consistency equa- 
tion (|^) must be considered as an operator equation. 
Thus, we are obliged to introduce the band bottom oper- 
ator u, which is a solution of Eq. rtq). A similar approach 
has already been used in Ref. p5j The charge operator 
Q becomes then a superoperator Q[u\. In the following 
we consider only the voltage range v a < v < v*, and in 
case of multiple stable states only the one correspond- 
ing to the charged well (denoted U3(v) in the preceding 
discussion) . 

Eq. (||) is non-linear. Instead of trying to find a gen- 
eral solution, we linearize it, writing u — u + 5u, where u 
is the solution of Eq. (g) for the average quantities. With 
the "parameter" of the expansion being a^a^ — (a^a^), 
we write the charge operator as follows, 



Q[u] = (Q[u]) + Qx + Q 2 , 

Qi = ou; Go u = — > 

e ou 



(20) 



Q 2 = e^^/ dcje- wt / dE z v af} {E z -u) 

x [al c (E z ,p ± )dp(E z + w,p±)- 
- (a 1 a {E z ,p ± )a l3 (E z + u,p±))] . 

Here the average charge iff) is given by Eq. (|l0|) and is 
a function of the average-potential it (not of the opera- 
tor u). The capacitances Cq[u] determines the charge 
increment in the well in response to an increment of 
the potential in the well. It is one of the important re- 
sponse functions which determines the new parameter A 
which eventually enters the Fano factor and was, there- 
fore, mentioned already in the introduction (see Eq. (1), 
U = u/e). On the stable branch (3 in Fig. 3) the "ca- 
pacitance" Co changes sign, and for the voltage v* it is 
negative and equal to — (Cl + Cr). The operator Q 2 
appears as the operator for the bare (unscreened) charge 
and is a function of the average potential u. For the 
operator Su we obtain 



Su(t) 



C L + C R + C [u] 



»(*)• 



(21) 



For v — » v* this expression diverges. We will show later 
that this divergence is at the origin of the noise enhance- 
ment as the voltage v* is approached. 



C. Current operator 

Following Ref. ^3[ we write the current density opera- 
tor in the contact a, 

Ut) = ^EE / dE z dE> z e« E *- E *»A% 7 (E z ,E>) 



0j p± 

xa\{E z ,p 1 _)a 1 {E' z ,px), 



(22) 



G 



with current matrix elements 

A% 7 (E Z) E' Z ) = S aP S ai - sl (E z )s^(E' z ) (23) 
which arc themselves functions of the operator u. Using 



Eq. (|17|) , we obtain 

A LL {E Z ,E Z ) = A RR {E Z ,E Z ) = -A RR (E Z , E z 

Y L [E Z )T R {E Z ) 



- -A LL (E Z ,E Z ) - {Ez „ Eo)2+T2/4 - (24) 

It is immediately seen that jji) = — (jr) and the average 
current (Jl) reproduces Eq. (|12|). 

The current operator ( p2| ) also must be linearized in Su 
in the same manner as we did it for the charge operator. 
We obtain 



ja[u] = (ja[u\) +jla +32a, 
djjg) , 

du 



e 



(25) 



32c 



2ttA 



EE/ dE z dE>e^ E '*- E ^A^(E z -E>) 



al(E z ,p±)a 1 (E' z ,p ± )- 



al(E z ,p 1 _)a 1 {E' z ,p ± ) 



The current operator j\ a is proportional to J ai the incre- 
ment in current at contact a in response to a potential 
variation in the well. J a is the second important response 
function which enters the Fano factor and has also been 
mentioned in the introduction (see Eq. (2)). j 2a is the 
bare current density and via the current matrix A 0J de- 
pends on u (not on u). Note that Jl — —Jr and that in 
the introductory section we have suppressed the contact 
index and have defined J = cAJl — —gAJr, U — u/e. 



D. Current-current fluctuations and the Fano factor 



Now we are in a position ta_evaluate the shot noise 
spectral power S a p, defined asEj 



2irS(u> + uj')S a i3{u)) = 



dtdt' e lujt+lut ' 



Ja(t)W) +h(t')j<x(t)) ' 2 (j a (t)) (5>(t' 



(26) 



We consider the zero-frequency limit and use the property 2 ^ 

(al(E 1 ,p±)ap(E 2 , p±)at f (E 3 ,p' ± )as{E 4: ,p'_ L )) - (d t a (E 1 ,p±)ap(E 2 ,p±)) (ci^Ea, p' ± )a s (E 4 , p' ± )) 
= S aS 8 fJj S p±pl 8(E 1 - E i )8{E 2 - B 3 )/«(^i +pl/2m) \l - f (E 2 + pl/2m) 



(27) 



After simple but tedious algebra where we replace at 
an intermediate stage ((E z - E ) 2 + T 2 {E Z )/A)- 2 by 
AttT-'\E )5{E z - E ), we obtain S LL = S RR = -S LR = 
—S R l (thus checking the conservation of current) with 



Sll = Sll(0) = 



2e 2 iy 2 T L (E )T R (E ) 



A r(E ) 3 

x (v - u + E F - Eo) {T 2 L (E ) + T 2 R (E ) 
eAJh 



+ 2 



Cl + C R + Co 

eAJ L 
Cl + C R + Co 



[T R (Eo)-T L (Eo)} 



(28) 



for v a < v < v* and Sll — otherwise. 

With the help of the spectral density of the cur- 
rent noise, evaluated on the stable branch of the I-V- 
characteristic, we can find the Fano factor. With the 
interaction energy 

eAJ L 



A 



Cl + C R + Cq 



Flm — Fl,r(Eo) 



we find 



F = 



2A [T R - T L ] + 2A 2 



[r L + r R ] 2 



(29) 



Eq. ( p9| ) is the main result of our work and together with 
the definitions of Jl and Co and with the equation (^|) for 
u(v) (the stable charged branch solution) in an implicit 
form gives the dependence of the noise power on the ap- 
plied voltage v. Using Eq. ( p9| ) for the Fano factor and 
the definition of the asymmetry Ar = (Tj, — r#)/2 gives 
the Fano factor, Eq. (||), mentioned in the introduction. 

The dependence on the Fano factor on the external 
voltage v is illustrated in an explicit manner in Fig. 5. 

There are several important points to comment on. 

First, the function A(v) accommodates the interac- 
tion parameter C e and represents the contribution of 
interactions to the noise power. For A = jjno in- 
teraction) we reproduce the free-electron resultu Fq = 
(r| + T R )(r L + T R )~ 2 , as already discussed in the in- 
troduction. Note, however, that due to the energy de- 
pendence of the tunneling rates for a symmetric dot, 
&l = o-Ri one has Tl < r#, which implies 1/2 < Fq < 1. 

Furthermore, for v — > v* the denominator C'l + Cr + 
Co(v) of the function A(v) quite generally diverges as 
(v* — v) -1 / 2 , while the function Jl(v) in the numera- 
tor stays finite. Consequently, the Fano factor diverges 
according to (v* — v) . This is clearly seen in Fig. 5. 



7 



For v = v a = C^{C L + C R ){E -E F ) the Fano factor 
can be calculated in a closed form. We have 



CoK) = c e 



which gives 



-ev 2 



Ll R 



r L + r 7 



v 



V V 



(d) 



C e T L T R 



{C L + C R )(T L + T R ) + C e T L 



< 0. 



It is seen that — T R < A(v a ) < 0, implying 1/2 < 
F(v a ) < 1. Thus, we have described the transition from 



F 




FIG. 5. The full Fano factor F (solid line) and its 
non-interacting part Fq (bold dashed line) as a function of 
applied voltage v. The parameters are chosen as follows: 
a L = a R ; C L = Cr; (a) C e = WC L , E = 3E F /2 (this set 
corresponds to that of Fig. 2-4); (b) C e = 10Cx, E = 3E F ; 
(c) C e = 4C L , Eo = 2E F ; (d) C e = C L , E = 3E F /2. The 
voltage v* (in units of Eo) is equal to 2.70 (a), 2.26 (b), 2.14 
(c), 2.02 (d). 

sub-Poissonian to super- Poissonian shot noise. The 
range of super-Poissonian noise (F(v) > 1) lies at volt- 
ages close to v*; this range is the larger the stronger the 
interaction is (the parameters C e /Cr, and C e /C R grow) 
and shrinks with increasing distance of the resonance 
from the equilibrium Fermi level, i. e. with increasing 
Ep I Eq. Note that the onset of the super-Poissonian noise 
(the voltage where F[v] — 1) is unrelated to the point of 
zero differential resistance (dj/dv = 0). Of course, these 
two points may lie quite close to each other. In addition, 
the quantity Jl changes sign from negative to positive 
when the voltage approaches v*. This implies that A 
also changes sign from negative to positive at the same 
voltage. 

Finally, with respect to the role of interactions, two 
scenarios emerge: To avoid cumbersome expressions, we 
perform the analysis for the symmetric case <zl = a R , 
Cl = Cr. Then for "weak interactions", 



C\ . C= If, I -1 



Ep 



we have F(v a ) < F (v a ). Thus, the curves F(v) and 
F a (v) cross each other only once, when A(v) = 0. This 
means that interactions suppress the shot noise for "low" 
voltages and enhance it for "high" voltages, eventually 
driving the noise power to super-Poissonian values for 
voltages close to v*. For "strong interactions", C e > C, 
we obtain F(v a ) > F (v a ), and the curves cross each 
other twice. In this case, interactions first enhance the 
shot noise, then suppress it, and then enhance it again. 
In the limiting case C e S> Cl one has F(v a ) = 1; then 
the shot noise starts from the Poissonian value. Fig. 5a 
and Fig. 5b represent the "strong interaction" scenario. 
Fig. 5d illustrates the case of "weak interactions" , while 
Fig. 5c shows the marginal case C e = C, where one has 
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F(v a ) = F (v a ). 



V. CONCLUSIONS 

We now discuss the assumptions we made in the course 
of the derivation of our result, Eq. (||). First, in order 
to calculate the shot noise, we linearized the charge and 
current operators in the fluctuations of the band bottom. 
Close to the voltage v*, at the linear stability thresh- 
old, where the solutions corresponding to the charged 
well vanish, these fluctuations become large, and the lin- 
earization is not justified any more. Since we evaluate 
only the zero- frequency response, the quantity u 2 (t) is 
not amenable to our treatment, and the limits of this 
approximation can not be determined. We emphasize, 
that in the crossover region from sub- to super-Poissonian 
noise our approach is well justified. Thus, our discussion 
can be used to determine the range of voltages in which 
super-Poissonian noise can be observed. 

To determine the effect of interactions on the Fano 
factor, we have neglected the finite width of the resonant 
level in the evaluation of A. For comparison with ex- 
periment the broadening of the resonant level (finite T) 
should be taken into account. This would smear all the 
curves; in particular, the noise power at v = v* would 
acquire a finite value. Moreover, if the tunneling rate T 
exceeds the hysteresis interval v* — Vf,, an extended range 
of differential resistance will be observed instead of the 
hysteretic behavior. We believe that is the case in the 
experiments of Ref. ||. Note also that in our treatment 
the interactions make the Fano factor discontinuous ai 
v = v a , where it jumps from the non- interacting valuetl 
Fo(v a ) to the renormalized one F(v a ). This discontinuity 
is also healed by the introduction of the finite width of 
the resonant level. 

In our work we have treated a single resonant level. 
It is clear that wells which permit to fill a number of 
subbands can be expected to exhibit even richer behav- 
ior than indicated in this work. Since interactions play 
a dominant role in the wide area resonant double bar- 
rier structure considered here, it might be interesting to 
analyze the noise power of a single discrete level in a 
quantum dot under high bias and ask if in fact a similar 
behavior results. It would also be interesting to investi- 
gate the effect of potential fluctuations not only in the 
well but also in the accumulation and depletion layers in 
the cathode and anode. In the terminology of Ref. |l6| we 
have in such a model internal degrees of the well which 
couple to the cathode. This would modify the current in- 
jection into the well (depending on the charge state) and 
possibly provides a mechanism which can lead to Fano 
factors smaller than F = 1/2. 

In conclusion, we investigated current-voltage charac- 
teristics and noise power in resonant tunneling quan- 
tum wells in the nonlinear regime. We treated the low- 
transparency case, when the total width T, due to tunnel- 



ing, is much lower than all other relevant energy scales. 
(For the evaluation of the interaction effects, we have 
taken the transmission coefficient as delta-function in 
energy). The interaction is taken into account via the 
charge accumulated in the well due to capacitive cou- 
pling to the reservoirs. In this regime, we discovered 
that interactions dramatically affect transport properties 
of the quantum well. They produce a range of volt- 
ages (vb < v < v*) where the well can be in one of 
two stable states, which correspond to a charged well 
and an uncharged well, respectively. For the current- 
voltage characteristics, this leads to a hysteretic behav- 
ior, as shown in Fig. 4. Interaction effects are even more 
pronounced in the noise power: At low voltages inter- 
actions can either suppress or enhance the noise power 
below or above the shot noise level predicted by a free- 
electron theory, and at large voltages the noise power 
diverges as (y* — v)^ 1 for v — > v*, and becomes super- 
Poissonian close to v* . Both the noise suppression below 
the non-interacting valueB'ta and the Strang enhancement 
of the noise to super-Poissonian valuestJU have been ob- 
served experimentally. We stress that the range of super- 
Poissonian noise is generally different from the negative 
differential resistance range, and may extend well outside 
the hysteresis range (see the example shown in Fig. 5d). 
We believe that our results provide an explanation of the 
role of interactions on the shot noise in resonant wells. 

Clearly, the interaction effects in the noise power which 
we have discussed are very interesting. While the search 
for interaction effects in mesoscopic physics represents a 
huge effort, there are up to now, at best very few clearly 
identified signatures for which both experiments and the- 
ory agree. A notable exception are of course the Coulomb 
blockade effects. We emphasize that the effects discussed 
here, occur in a state which is a far from equilibrium 
transport state. Resonant tunneling systems have the 
important advantage that they represent, at least con- 
ceptually, a very simple system and thus provide a unique 
testing ground for both theories and experiments. 
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